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\q ' We present two-dimensional spectroscopy covering the rest-frame wavelengths of 

, strong optical emission lines in six luminous submillimetre galaxies at z = 1.3-2.5. 

Using this near-infrared integral field spectroscopy together with Hubble Space Tele- 
scope ACS and NICMOS imaging we map the dynamics and morphologies of these 
systems on scales from 4-11 kpc. Four of the systems show multiple components in 
their spatially-resolved spectra with average velocity offsets of ~180kms _1 across 
\ 8 kpc in projection. From the ensemble properties of eight galaxies, from our survey 

c/3 . and the literature, we estimate the typical dynamical masses of bright submillimetre 

galaxies as 5±3xlO n M Q . This is similar to recent estimates of their stellar masses 
- suggesting that the dynamics of the central regions of these galaxies are baryon 
dominated, with a substantial fraction of those baryons in stars by the epoch of ob- 
' servation. Combining our dynamical mass estimates with stellar luminosities for this 

^ . population we investigate whether submillimetre galaxies can evolve onto the Faber- 

Jackson relation for local ellipticals. Adopting a typical lifetime of Tburst ~ 300 Myr 
for the submillimetre-luminous phase - using the latest estimates of gas masses, star 
formation rates and AGN contribution to the bolometric luminosities - we find that 
the stellar populations of submillimetre galaxies should fade to place them on the 
Faber-Jackson relation, at Mr- ~ —25.1. Furthermore, using the same starburst life- 
time we correct the observed space density of submillimetre galaxies for the duty cycle 
to derive a volume density of the progenitors of ~ 1 x 10~ 4 Mpc~ 3 . This is consistent 
with the space density of local luminous early-type galaxies with Mk ~ —25.1, indi- 
cating that submillimetre galaxies can evolve onto the scaling relations observed for 
local early- type galaxies, and the observed population at z ~ 2 is then sufficient to 
account for the formation of the whole population of > 2>L* K ellipticals seen at z ~ 0. 

Key words: galaxies: high-redshift, submillimetre; — galaxies: evolution; — galaxies: 
star formation rates, AGN; — galaxies: individual 



1 INTRODUCTION 

Deep optical and near-infrared imaging with the superla- 
tive resolution of the NICMOS and ACS cameras on- 
board Hubble Space Telescope (HST) has made it pos- 
sible to study the morphologies and colours of high- 
redshift far-infrared luminous galaxies identified in the 



sub-millimetre (submm 


waveband in unprecedented lev- 


els of detail l|Smail et al. 


|l998 


l200llChaDman et alJl2003b: 


ISmail et all2004t IPone t 


t alJI 


200,4 lAlmaini et all200ffl. By 



combi ning this data with spectrosco pically confirmed red- 
shifts dChapman et al.ll2003al . l2005f) . we are beginning to 
understand the processes which triggers the immense bolo- 
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metric luminosity output in these galaxies, allowing us to 
address issues such as their true contribut ion to the star- 
forma tion rate history of the Universe ('e.g.. lChapman et alJ 
120051) . Submm galaxies (SMGs) are amongst the most 
bolometrically luminous galaxies in the Universe (with 
Lboi ~10 13 Lq and a median redshift of < z >~2.2). 
However, in contrast to quasars at the same early epochs, 
their bolometric output appears to be dominated by star- 
formation rather than AGN activity ([Swinbank et aL 2004; 
I Alexander et alJl2005|) . It is clearly important to establish 
firmly the masses for these galaxies in order to understand- 
ing the rapid evolution of the submm population. Such diag- 
nostics will allow us to determine how they relate to present- 
day galaxies and in particular to test whether they represent 
the formation phase of the most ma ssive elliptical ga laxies 
at the present-day, as many suspect jLillv et al] ll999l. 

To probe the dynamical structures of these frequently 
complex systems requires a reliable separation of the spa- 
tial and spectral information. One particularly powerful 
approach to achieve this is using millimetric CO emis- 
sion line maps to trace the distribu tion and kinemat- 
ics of dense gas within the galaxies Fraver et al.l Il999l 



20031: iNeri et alJl2003l: iGreve et alj|2005l) especially a t high 



resolution (lOenzel et al . 2003; Downes fc Solomon! |2 
iTacconi et aljl200fj) . Unfortunately, prior to the commis- 
sioning of the Atacama Large Millimetre Array, this ap- 
proach remains observationally very demanding, typically 
requiring 30 hours per source with the IRAM interferome- 
ter and can only be applied to the most massive and gas-rich 
SMGs. An alternative tool which is less demanding in terms 
of telescope time and can provide spatially-resolved dynam- 
ics of SMGs with sufficient spatial and velocity resolution to 
study their internal structure is near-infrared integral field 
spectroscopy targeting the rest-frame op tical emission lines 
such as [Om1A5007 a nd HaA6563 (e.g. iTecza et al]l2004t 
ISwinbank et~ai]|2005l) . The spectroscopic maps so produced 
allow us trace the dynamical and structural properties of 
SMGs on scales of a few kpc, pin-point the sites of active 
star formation and identify non-thermal emission from ac- 
tive galactic nuclei (AGN) in components within these sys- 
tems. However, they are also subject to potential biases due 
to extreme extinction within some regions of the SMGs and 
from o utflows in the emission-line gas which is being ob- 
served dBower et all 12004 IWilman et alll200Sh . While the 
star formation rates estimated from the Ha are typically 
a factor of 10-100 x less than that predicted from the far- 
infrared emis sion (suggesting t o 1-4 magnitudes of dust ab- 
sorption; q.v. lSmail et al.l2004h . it is clear that the Ha emis- 
sion provides a more secure probe of the star-formation rates 
and dynamics than either the Lya emission or rest-frame 
UV continuum and can thus provide a reliable and observa- 
tionally efficient probe of the activity and dynamics of this 
important population. 

Using the UIST near-infrared Integral Field Unit (IFU) 
at UKIRT and the GNIRS IFU at Gemini-South we have 
studied the rest-frame optical emission line properties of six 
submillimetre galaxies at z = 1.3-2.5. These observations, 
together with high-resolution imaging from HST allow us 
to probe the physical cause of the rapid evolution of sub- 
millimetre galaxies, enabling us to test whether these far- 
infrared luminous galaxies comprise merging systems which 
are likely to be the progenitors of local massive ellipticals, or 



whether instead they are simply high-luminosity episodes in 
the history of more mundane galaxies. Our paper is laid out 
as follows: in f|21 we present the observations and their re- 
duction. In S^we discuss the analysis of these data, while in 
S0|and we present our discussion and conclusions respec- 
tively. We use a cosmology with Ho = 70kms~ 1 , Qm = 0.3 
and Q,a ~ 0.7 in which 1" corresponds to 8.2 kpc at 2 = 2.5 
and 8.5 kpc at z — 1.5. 



2 OBSERVATIONS AND REDUCTION 

Th e SMGs in this pap e r come from surveys at 85 or 1200^tm 
bvlScott et alJ J2002D : IGreve et all (l2004h and IWebb etai] 
J2003bl). These s ources were s pectroscopically identified by 
IChapman et aL <2003al |2005|) a nd co nfirmed with precise 
redshifts from Swinbank et alJ 12 0041 based on accurate 



radio positions (e.g. llvjsmi^t^dTl2002T) . We also include 
the recently discovered hyper-luminous far-infrared-s e lected 
galaxy MIPS J142824.0+352619 from iBorvs et al l (12006ft 
aswell as an optically faint radio galaxy (OFRG) which was 
tentatively d etecte d as a millimeter continuum source in 
IGreve et al.l i2005T) for which a spectroscopic redshift was 
measured by Chapman et al. (2006). The coordinates and 
redshifts for each source are reported in Table 1 and we 
briefly review the properties of each individual source here. 

2.1 Sample 

CFRS 03.15: Fir st identified in th e Canada-UK Deep 
Submm Survey bv lWebb etaD <l2003aft . this submm galaxy 
has 5850 = 4.4 ± 1.1 mjy and a 1.4GHz flux density of 
<5i.4GHz = 226 ± 12/jJ y. Sp ectroscopic obse r vation s from 
IChapman et alJ i2005l) and ISwinbank et alJ (I2004D yield 
a redshift of z = 1.4076 which implies a luminosity of 
LF/fl =5.8±1.5x 10 12 L p. Th e rest-frame UV spectroscopy 
from IChapman et al] i2005T) shows strong C IV emission, 
which along with the high [Nil] /Ha line ratio and apparent 
underlyin g broad Ha line in the near-infrared spectra from 
Swinban k et al] B2004f) suggests substantial AGN activity. 

N2 850.7: This submm-selected galaxy has 
5 , 85o=9.0±2.4 mjy, lies at 2 = 1.488, and h as a far-inf rared 
lumin osity of Lfir = 6.4 ± 1.7 x 10 12 L© llChapman et al] 
2005). The existing optical and near-infrared spectroscopy 
of this source suggest that the bolometric luminosity is 
predominantly arising from a star-b urst, with low [Nil] /Ha 
and Civ/Lyq emissi on line ratios dChapman et al.l 120051 : 
ISwinbank et al] l2004h . Interferometric observations of the 
molecular CO emission in this system bv lGreve et al] ]2005h 
failed to detect the CO(2-l) line. If the gas reservoir within 
this galaxy is close to the systemic redshift, indicated by 
the Ha emission line, then this non-detection implies an 
upper limit of < 1.8 x 10 10 Mq on the H2 gas mass. 

N2 1200.18: The near-infrared s pectroscopy of 
N2 12 00.18 indicates a redshift of z = 2.592 ISwinbank et al] 
l2004h . and its 1200/im flux density of 2.2±0.6 pJy 
then corresponds to a far-infrared luminosity of 
10.9 ± 3.0 x 10 12 L for a reasonable tem- 
plate SED. The radio flux density (Si.4GH* = 180 ± 22/iJy) 
is relatively strong, suggesting a possible contribution 
from an AGN. This is confirmed by the existing longslit 
optical and near-infrared spectroscopy which shows high 
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Figure 1. Left: True colour HST ACS/NICMOS B450-/gi4iii6o-baiid image of CFRS 03.15 at z — 1.4076 showing the complex morphology (the colour image is displayed 
in log scale to emphasise the morphology). Middle: HST NICMOS -ffigp-band image of CFRS 03.15 with the Ho emission line map from the ONIRS IFU observations 
overlaid as contours, showing that the bright knot located ~1.3 (~11 kpc) to the South- West is associated with the system. The solid bar represents the FWHM of the 
seeing disk for the IFU observations. Right: Spectra from the two components seen in the HST imaging of CFRS 03.15 from the Kcck/NIRSPEC observation. The top 
spectrum is offset in flux for clarity. From our spectroscopy, the compact component to the South-West is blue-shifted by 90±20kms — 1 in projection from the nucleus, 
suggesting they arc two components whose interaction has promoted the strong starburst and disturbed morphology. It is interesting to note that both spectra have 
strong [Nlll]/Ha emission line flux ratios which may indicate that both components host active AGN. 
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Figure 2. Left: A true colour image of N2 850.7 comprising HST ACS igi4-band and NICMOS iJigp-band images illustrating the complex morphology of this system. 
Middle: NICMOS -ffigQ-band image with Ha emission line intensity overlaid as contours. Each panel is 4 arcscc square (34 kpc at z — 1.49) and has East left and North 
top. The bar shown on the NICMOS image indicates the 0.6" resolution of the IFU observations. The system comprises at least three components, marked as A, B and C 
(and there is also tentative evidence for weak Ha emission from the source ^ 1 to the West, suggesting it may be associated with the system). The continuum colours of 
C are bluer than the other components. The intensity of the Ha emission traces the J/^eo-band morphology and it is interesting to note that the star formation appears 
distributed across the SMG. Right: Spectra covering the rcdshiftcd Ha emission in the three components of N2 850.7 from the UIST IFU observations. The spectra have 
been offset in flux scale for clarity. The lower panel shows a sky spectrum. By ccntroiding the emission lines, we find velocity offsets between components A & B and A 
& C of 215 zb 80 kms -1 and 135 ± 90kms _1 in projection respectively. 
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Figure 3. Left: HST igi4-band image of N2 1200.18 with the [Olll] A5007 and Ha emission line intensity maps from the UIST IFU observations overlaid as contours. Both 
panels arc displayed on a log scale to emphasis the morphology. The HST image shows a a central bright, compact component and a lower-surface brightness extension to 
the North. The [OlIl]A5007 and Ha emission also have different morphologies, with the strong, compact [OlIl]A5007 probably arising from (unresolved) AGN activity in 
the bright central component, whilst the Ha morphology appears resolved and extends across the whole system. Right: Spectra from the central component and northern 
extension of N2 1200.18 around the [OlIl]A5007 and Ha emission lines from the UIST IFU observations. The spectrum of the central source has an [OlIl]/Ha emission line 
flux ratio of 2.0i0.5 which, combined with the unresolved [Olll] emission suggests AGN activity. The northern component has a velocity offset of 250±80kms — from 
the nucleus. This northern component has star-burst (rather than AGN) characteristics, with a lower [OlIl]/Ha and [Nll]/Ha emission line flux ratios. Bottom: The OH 
airelow emission from the skv (continuous line), as well as the Mauna Kca atmospheric transmission (dotted line) arc also shown in the lower panel. 
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Figure 4. Left: HST Igi4-band image of N2 1200.90. The galaxy has a compact morphology (-g^O.2 (1.6 kpc) FWHM) in this continuum imaging. Right: Spectrum of 
N2 1200.90 from the UIST IFU observations showing Ha at a redshift of z — 2.198. The rest-frame UV spectrum from Chapman et al. (2006) shows clear signs of AGN 
activity, yet the apparently low [Nil] /Hex emission line ratio from our near-infrared spectroscopy is more suggestive of star-formation. However, we note that the Ha 
emission line width of 500+100 km s —1 may indicate (non-thermal) AGN activity. 
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Figure 5. Images of MIPS J142824. + 352619 generated from the UIST IFU observations. Left: JC-band continuum image (collapsed from the datacube between 2.0 and 
2.4/xm) with contours overlaid from the Jf-band (collapsed between 1.5 to 1.8/im). Middle: Continuum subtracted Ha image. Right: X-band image with contours from 
the continuum subtracted, Ha narrow-band image overlaid. The "X" in each panel marks the position of the center of the .ff-band light distribution. In the final panel 
we show the collapsed, one-dimensional spectrum around the redshifted Ha emission from MIPS J142824. 0+352619. The Ha redshift of z — 1.328 + 0.003 is in excellent 
agreement with the CO(5-4) emission, and the Ha line width and emission line flux are also comparable to those given in Borys et al. (2006). 
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Figure 6. Left: HST Jg 14 -band image of SA 22.96 with the Ha emission from the UIST IFU observations overlaid as contours. Right: Collapsed spectra around the Ha 
emission from the Western component (A) from the UIST IFU observations of SSA22.96. We identify A as the submm source at redshift z — 2.517 through identification 
of strong Ha and [Nil] emission at the same redshift as derived from the rest-frame UV spectroscopy in Chapman et al. (2005). We note that there is tentative evidence 
for weak Ha emission from the structure ~ 1.5" to the East, suggesting it may be associated with the system. 
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Table 1. 
Log of IFU Observations 



SMG 




RA 


Dec 


2 


texp 


UIST/UKIRT 












SMMJ163639.01- 


f 405635.9 (N2 850.7) 


16:36:39.40 


+40:56:38.0 


1.4880 


21.6 


SMM J163655.80- 


f 405914.0 (N2 1200.18) 


16:36:55.90 


+40:59:12.0 


2.5918 


21.6 


RG J163655. 05+410432.0 (N2 1200.90) 


16:36:55.05 


+41:04:32.0 


2.1981 


14.4 


MIPS J142824.04 


-352619.0 


14:28:24.40 


+35:26:21.8 


1.328 


21.6 


SMMJ221804.42- 


f 002154.4 (SSA 22.96) 


22:18:04.43 


+00:21:53.2 


2.517 


10.8 


Gemini/GNIRS 












SMM J030227.73- 


f 000653.5 (CFRS 03.15) 


03:02:27.6 


+00:06:52.5 


1.4076 


7.2 



Table 1. Note: The submm names/position give the 850/jm centroid whilst the RA and DEC marked denote the 1.4GHz radio centroid. 



[Nil] /Ha and Civ/Lya emission line ratios as well as a 
broad component underlying the Ha emission line (with 
FWH M of 2900 ± 400 km s" 1 ) from a partially obscured 
AGN iChapman et alJl2005t ISwinbank et alJl2004D . 

N2 1200.90: The 1.4-GHz map of the ELAIS-N2 field 
shows a slightly resolved source around this galaxy, with 
a peak flux density of 47 ± 10/iJy beam" 1 and a total flux 
density of 63±21/iJy (Biggs & Ivison 2006). A spectroscopic 
redshift z — 2.195 was measured through the identification 
of Lya and C iv emission lines, the ratio of these suggests 
weak AGN activity is possible (Chapman et al. 2006). 

MIPS J142824. 0+352619: Far-infrared imaging 
with the Spitzer Space Telescope recently unco vered an ap- 
paren tly hyper-luminous galaxy at z ~ 1.325 feorvs et alJ 
2006). Multi-wavelength follow-up of this galaxy at 
submm wavelengths suggest a far-infrared luminosity of 
3.2±0.7xl0 13 Lg. Unlike most other z > 1 sources of 
comparable luminosity, MIPS J142824. 0+352619 lacks any 
trace of AGN activity, and therefore the immense lumi- 
nosity appears to arise due to a ^SOOOMQyr" 1 star-burst. 
However, spectroscopic observations in the z-band have 
shown that this galaxy is aligned with a foreground 0=1.03 
elliptical galaxy, suggesting the possibility of gravitational 
amplification of the background source. This would provide 
a much more mundane explanation for the apparently 
immense luminosity of this galaxy. 

SSA 22.96: Rest - frame UV spectroscopy of SSA 22.96 
from Cha pman et alJ Q2Q05) indicates a redshift of z = 
2.517 from UV absorption features, with little indica- 
tion of AGN activity. This radio pre-selected galaxy has 
850 fim and 1.4 GHz flux densities of S85o=9.0±2.3mJy and 
Si.4=43.8±10.1/xJy which suggests a far-infrared luminosity 
of Lfir ~8xl0 12 L . 

Deconvolving the exact contributions from the star for- 
mation and AGN activity within these gal axies is di fficult 
using just the UV slit spectroscopy from IChapman et alJ 
(2005). Ideally we need to use the well-developed spec- 
tral indicators based on res t-frame optical emission lines 
JVeilleux fc Osterbrocklll987l) which fall in the near-infrared 
for these high-redshift galaxies dSwinbank et al.l l2004). Go- 
ing one step further, and coupling the spatial coverage from 
an IFU with coverage around the rest-frame optical emission 
lines allows us to locate and isolate the components hosting 
the AGN in these systems, determine dynamical masses, as 
well as searching for extended halos and/or companions. 



2.2 UIST Near-Infrared Integral Field 
Spectroscopy 

Spectro-imaging observations of five targets were made 
with the UIST near-infrared integral field spectrograph on 
UKIRT. The UIST IFU uses an image sheer to take a 
3.3"x6.0" field and divides it into 14 slices of width 0.24". 
The dispersed spectra from the slices are reformatted on the 
detector to provide 2-D spectro-imaging, in our case using 
the HK grism, at a spectral resolution of A/AA ~ 1000 
and covering a wavelength range of 1.4-2.4 /im. N2 850.7 
and N2 1200.18 we observed between 2004 August 20 and 
September 04, whilst the remaining targets were observed 
between 2005 June 25 and June 30 . All observations we 
taken < 0.6" seeing and photometric conditions. Observa- 
tions were carried out in the standard ABBA configuration 
in which we chopped by 12" to blank sky to achieve good 
sky subtraction. Individual exposures were 240 s and each 
observing block was 7.2 ks which was typically repeated two 
or three times, thus the total integration time for each object 
were 10-20 ks (Table 1). 

To r educe the data we u sed the relevant orac-dr 
pipeline llCavanagh et alJ 120031) which sky-subtracts, ex- 
tracts, wavelength calibrates, flat-fields, flux calibrates and 
forms the datacube. To accurately align and mosaic the in- 
dividual datacubes together we created white light (wave- 
length collapsed) images around the redshifted Ha emission 
line from each observing block and used the peak intensity to 
centroid the object in the IFU datacube. We then spatially 
aligned and co-added the individual data-cubes to create the 
final mosaic. 



2.3 GNIRS Integral Field Spectroscopy 

We used the GNIRS IFU spectrograph to observe the z = 
1.4076 submm galaxy CFRS03.15 2 . The GNIRS IFU uses 
an image sheer to take a 3.2" x 4. 8" field and divides it into 
21 slices of width 0.15". By re-arranging the slices end-to- 
end at the detector, the IFU allows images to be dispersed 
into full length spectra whilst preserving the complete 2-D 

1 Programme ID: U/04A/162 & U/05A/45: The United Kingdom 
Infrared Telescope is operated by the Joint Astronomy Center 
on behalf on the UK Particle Physics and Astronomy Research 
Council. 
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spatial information content. This allows 3-D spectroscopy 
across a contiguous 3.2" x4. 8" field at 0.f5" per pixel. We 
used the //-band filter with the 32 lines/mm grating which 
results in a spectral resolution of A/AA ~1700 and records 
a wavelength coverage of between 1.37 and 1.87/xm. Obser- 
vations were taken on 2005 October 20 in photometric con- 
ditions, but moderate (~1.0") seeing. We used the same 
observing strategy as with UKIRT: observations were car- 
ried out in the standard ABBA configuration in which we 
chopped away to sky by 12" to achieve good sky subtrac- 
tion. Individual exposures were 600 s seconds and each of the 
ABBA observing blocks was repeated three times to give a 
total integration time of 7.2 ks. To flux calibrate the data, 
we observed the UKIRT standard star FS10. We used the 
Gemini IRAF package to extract, flat-field and wavelength 
calibrate the data and used idl to build and combine the 
datacubes. 



2.4 HST Optical and Near-Infrared Imaging 

HST Advanced Camera for Surveys (ACS) observations of 
N2 850.7, N2 1200.18 were obtained from the HST pub- 
lic archive 3 (Program ID #9761) whilst ACS images of 
SSA 22.96 was observed as part of Program #10145. The 
data consist of dithered exposures with the F814W fil- 
ter, taken in lowsky conditions using the default four- 
point ACS-WFC-dither-box configuration. This pattern en- 
sures optimal half-pixel sampling along both coordinates. 
The typical integration time were 4.8 ks. We reduced the 
data using the latest ve rsion of the multidrizzle software 
flKoekemoer et al"ll200-f) using the default parameters with 
pixfrac=1 and SCALE=1. The resulting images have 0.05" 
pixels and are free from artifacts. Further details of these 
observations are gi ven in Borys et al. (2006b) (see also 
lAlmaini et"al"ll2005l) . 

HST WFPC2 observations of CFRS 03.15 were also ob- 
tained from the HST public archive. The data were taken in 
the F814W and F450W filters for total of 1.4 ks each (Pro- 
gram ID #05996 and #06556). These observations were re- 
duced using the standard STSDAS package in iraf. 

HST NICMOS observations of N2 850.7 and 
CFRS 03.15 were obtained in Cycle 12, and the tar- 
gets were observed using the NIC2 camera in the F160W 
filter for a total of 2.3 ks (Program ID #9856). We 
employed the standard four point spiral dither pattern, 
lowsky conditions and used the multiaccum read-mode. 

2 Programme ID: GN-2005B-Q-60: Based on observations ob- 
tained at the Gemini Observatory, which is operated by the As- 
sociation of Universities for Research in Astronomy, Inc., un- 
der a cooperative agreement with the NSF on behalf of the 
Gemini partnership: the National Science Foundation (United 
States), the Particle Physics and Astronomy Research Coun- 
cil (United Kingdom), the National Research Council (Canada), 
CONICYT (Chile), the Australian Research Council (Australia), 
CNPq (Brazil) and CONICET (Argentina) 

3 Obtained from the Multimission Archive at the Space Telescope 
Science Institute (MAST). STScI is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA con- 
tract NAS5-26555. Support for MAST for non-HST data is pro- 
vided by the NASA Office of Space Science via grant NAG5-7584 
and by other grants and contracts. 



Each exposure was corrected for a pedestal offset, and then 
mosaiced using the CALNICB task in iraf. Unfortunately the 
observation was affected by the South Atlantic Anomaly 
(SAA), and extra processing steps were required 4 . The 
final images appear very flat and have very low cosmic ray 
contamination. The reduction and analysis of these data is 
described in more detail in Borys et al. (2006b). 



3 ANALYSIS & RESULTS 

3.1 Velocity Determinations 

To accurately determine the redshift (or velocity offsets be- 
tween components) in the spectra, we identify and fit both 
the continuum level and emission lines (e.g. Ho?A6562.8, 
and [Nll]AA6548. 1,6583.0 where appropriate) simultaneously 
with a flat continuum plus Gaussian profiles using a y 2 fit 
and taking into account the increased noise in regions of 
strong sky emission. The error bars on the emission-line 
wavelength in each spectrum are derived by allowing the 
wavelength of the best-fit Gaussian profile to vary and allow- 
ing the signal to drop by A\ 2 = 1. The quoted error-bars are 
therefore conservative and reflect the signal-to-noise of the 
data. The velocity error is propagated through to the mass 
determinations in ^\ However, since the signal-to-noise in 
some of our observations is modest we also attempt to fur- 
ther quantify the error-bars on the velocity offsets by cre- 
ating artificial emission lines at the resolution of the UIST 
spectrograph. We produce spectra with emission lines with 
an intrinsic width of 150 kms -1 and add Gaussian random 
noise. We then refit the Gaussian profiles to the data, mea- 
suring the velocity offset betwen the input and output spec- 
trum. Using 10 4 simulations at each signal-to-noise, the rms 
velocity error between the measured fit and the input ve- 
locity as a function of signal-to-noise (S/N) is: 140 kms -1 
at S/N=3, 100 kms" 1 at S/N=5, 65 kms" 1 at S/N=7 and 
50 kms" 1 at S/N=9. As expected, these estimates are com- 
parable to the quoted velocity errors derived using our \ 2 
estimate. 

Before we discuss what can be learnt from our observa- 
tions as a whole, we briefly review the results from each of 
the six galaxies. 

3.2 CFRS 03.15 

In Fig. 1 we show a true colour 73/77-band HST image of 
CFRS 03.15. The galaxy shows a compact source surrounded 
by a diffuse halo of material in a ring-like structure. The 
colours of the nucleus (-B450 — 7814 = 1.8 7si4 — 7/160=2.5) are 
consistent with a star forming galaxy at z — 1.407, although 
the "ring" of material surrounding the nucleus is much red- 
der in (75 - I) Z 3.5, but bluer in (7-7/) = 1.8 ± 0.2, 
suggesting that [Oll]A3727 emission (A bs=8973A) may be 
contributing significantly to the /-band flux. There is also a 
compact knot (or companion) approximately 1.5" (12.5 kpc) 
to the South- West with similar colours to the nucleus. 

From our datacube we construct a wavelength-collapsed 
Ha image (collapsed between —300 and +300 kms" 1 from 

For a full description, see 

://www. stsci.edu/hst/nicmos/tools/post_SAA_tools. html 



http 
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the systemic redshift of IChapman et aDl2005l) and overlay 
this on the NICMOS H-band image (Fig. 1). The H-band 
and Ha morphologies are well matched, and we confirm that 
the bright knot located ~1.3" (~llkpc) to the South- West 
is associated with the galaxy (the extended Ha morphol- 
ogy and Ha velocity shear was al so tentatively detected 
in the Keck/NIRSPEC spectrum in lSwinbank et all J2004h 
in which the longslit was placed across both components). 
Since the observations were taken in ~1" seeing, we can- 
not discuss the small scale properties of each component in 
detail (or confirm the nature of the low surface brightness 
halo of diffuse material surrounding component A), however, 
these will be discussed in Blain et al (in prep) who use Laser 
Guided Adaptive Optics assisted i/-band Keck/OSIRIS IFU 
observations of this system. 

By extracting spectra from around the two brightest 
components we identify a velocity offset between the nu- 
cleus and the knot of 90±20kms~ 1 . We note that the 
[Nil] /Ha emission line flux ratios of components A and B 
are both high (0.95±0.05 and 1.05±0.05 respectively), sug- 
gesting that both components may host active AGN. In or- 
der to further search for signatures of non-thermal activity, 
we attempt to fit each of the spectra with an underlying 
broad-line component, but the improvement in % 2 over just 
fitting [Nll]+Ha+[Nll] is not significant. 



3.3 N2 850.7 

We detect Ha emission at 1.6 30/im (in agreement with 
the UV and Ha redshift from IChapman et alJ 120051 and 
ISwinbank et alJ l2004) and construct an Ha image of 
N2 850.7 by collapsing the datacube between 1.620 and 
1.640jitm. In Fig. 2 we show the HST NICMOS image of 
N2 850.7 and overlay on this the Ha image from the IFU 
observations. The Ha morphology of the central compo- 
nents is spatially extended and therefore to search for ve- 
locity sub-structure, we extract spectra from the regions in 
the IFU datacube around the brightest components seen 
in the HST imaging. From each individual component, the 
signal-to-noise of the Ha emission is only modest; never- 
theless, it is sufficient to centroid the Ha emission line in 
all cases and so search for velocity offsets between compo- 
nents. We identify a velocity offset between components A 
and B of 215 ± SOkms" 1 and an offset between A and C 
of 135 ± 90kms-\ The HST imaging and velocity offsets 
suggest that at least two galaxies are undergoing an inter- 
action, although whether component C belongs to the same 
structure as B (or A) will have to wait for higher resolution 
(adaptive optics) IFU observations. We note t hat the con- 
tinuum colours of A, B and C are very similar jSmail et alJ 
12004) . 

Using the Ha luminosity as a star-formation rate indi- 
cator, we derive star-formation rates of 90±20, 70±20 and 
60±20 M^yr -1 from components A, B and C respectively 
llKennicuttlll998h . In comparison, the star- formation rate 
for the whole system from the far-infrared luminosity is 
1100 ± SOOMfflvr -1 , im plying ~ 2 magnitudes of extinction 
(q.v. ISmail et al.ll2004h . We reiterate that neither the emis- 
sion line flux ratios from the I FU observations nor th e rest- 
frame UV spectroscopy from IChapman et all J2005D show 
any signs of AGN signatures. 



3.4 N2 1200.18 

The HST / 8 i4-band ACS image of N2 1200.18 (z = 
2.5918) appears complex with a bright knot and a low- 
surface-brightness extension or companion approximately 
0.7" (6kpc) to the North. Since the IFU observations of 
this galaxy cover H/3, [Oni]A5007 and Ha, we construct 
wavelength-collapsed images around the [Oni]A5007 and Ha 
emission lines and overlay these on the NICMOS -ffi60- 
band image in Fig. 3. The strongest [Om]A5007 emission 
appears to come from an unresolved point source in the cen- 
tral, bright 7gi4-band continuum component. In contrast, 
the Ha is much more extended and has features which 
match the extended Northern component seen in the 7si4- 
band. This suggests that the compact continuum source may 
host an AGN and indeed, both the strong, spatially unre- 
solved [Oni]A5007 and high [Om]A5007/Ha(=2.0±0.5) and 
[Nll]/Ha(=0.7±0.1) emission line ratios suggest AGN activ- 
ity. We note that the Ha/H/3 emission line flux ratio of the 
central and northern components are 8±3 and >5 respec- 
tively. 

We can use the emission line flux ratios of [OlIl]/H/3 
and [Nn]/Ha of N2 1200.18 to classify the two components 
seen in the UIST IFU data. The northern component has 
an [Nll]/Ha emission line flux ratio of 0.6±0.1, and an 
[Om]A5007/H/3 ratio of 0.5±0.2, together these place this 
compone nt it in the LINER r egion of the diagnostic dia- 
grams of iBaldwin et all dl98lTl . The high [Om]A5007/Ha 
and [Nil] /Ha emission line flux ratios from the central un- 
resolved component, indicate a Seyfert-2 type AGN. 

Using the velocity derived from the Ha emission, we find 
that the unresolved central component is redshifted from 
the northern extended component with a velocity offset of 
250 ± 100 kms -1 across 0.8" (~7kpc) in projection. 

Using deep Chandra observations of the ELAIS N2 
field it is also possible to compare the X-ray and spectral 
properties of N2 1200.18. This galaxy is detected in the 
hard (2-8 keV) X-r ay image with a flu x of 14.7 ± 2.1 x 
10" 15 ergs" 1 cm -2 jManners et al.ll2003ll We convert the 
observed 2-8 keV flux t o a rest frame 2-10 keV luminosity 
using the equations from I Alexander et all J2003bfl assuming 
a spectral index F=2. We derive Lx=9.1±1.3x 10 44 ergs _1 . 
Under the assumption that the [Om]A5007 emission 
line and the hard X-ray fluxes are isotropic, they can 
be used to investigate the intrinsic power of Seyfer t 
galaxies dMulchaev et al.lll994l ; lAlonso-Herrero et al.11199^1 . 
N2 1200.18 sits comfortably in the scatter of t h e \j\n jin- 
hx plot for local Seyferts from iMulchaev et alJ Jl994f) . al- 
though it lies at the high luminosity end, furthering the 
suggestion that N2 1200.18 hosts a high luminosity Seyfert 
2 AGN, and also suggesting that increased AGN activity 
results in increased photo-ionization of the surrounding ma- 
terial. This source, along with low-resolution spectroscopy 
of the [Olll]5007 emission from a larger sample of submm 
galaxies, is discussed in detail by Takata et al. (2006). 



3.5 N2 1200.90 

The HST lgi4 image of this galaxy shows a compact mor- 
phology (£0.2" (1.6 kpc) FWHM), which is matched by the 
Ha morphology from the UIST IFU observations. The Ha 
emission from this galaxy is spatially unresolved and we de- 
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tect no significant velocity structure across the galaxy. The 
line width of the Ha emission is 420±100 kms -1 and the in- 
tegrated line flux of 0.5x lO" 1 ^ Wm~ 2 then suggests a SFR 
of 15O±5OM yr- 1 llKennicut jll99l) . 



3.6 MIPS J142824. 0+352619 

The discovery that there is a fortuitous alignment of a mas- 
sive foreground elliptical galaxy at z=1.03, s uperimposed 
onto the z= 1.325 far- infrared- luminous galaxy feorvs et alJ 
l2006f) suggests the possibility that MIPS J142824. 0+352619 
may be gravitationally lensed. To this end, we tar- 
geted the Ha emission with the UIST IFU in order 
to search for morphological signatures of strong lensing 
which might explain the apparent immense luminosity of 
MIPS J142824.0+352619. 

In Fig. 5 we show the continuum subtracted Ha emis- 
sion line map of MIPS J142824.0+352619 overlaid on the 
K-band continuum image generated from the IFU. The Ha 
emission is spatially extended, across ~1" in projection, 
however, the H and Jf-band light is much more widely dis- 
tributed, showing likely contamination by light from the for- 
ground lens. The Ha emission appears elongated, but with 
an alignment which is not naturally ascribed to gravitational 
lensing (we would expect the brightest images to straddle 
the long-axis of the foreground lens potential). To search for 
velocity structure, we extract a series of independent spec- 
tra from the datacube along the arc, and place a limit of 
£125 kms" 1 on any possible velocity shear. We show the 
collapsed, one-dimensional spectrum around the Ha emis- 
sion in Fig. 5. The lack of velocity structure may be a further 
indication that the Ha morphology arises from two merging 
images of the same background source, with the third (much 
lower surface brightness) image undetectable to the South 
West. Whilst such a lens model is somewhat contrived, it 
would potentially yield an amplification ^10 for the back- 
ground source. 



3.7 SSA 22.96 

The HST ACS J-band imaging of SSA 22.96 shows several 
components around the submm galaxy (labelled A). Com- 
ponent A is the radio counterpart in the 1.4-GHz map from 
IChapman et alJ i2004f) . To align the IFU datacube with the 
optical imaging we construct a Ha image of SSA 22.96 by 
collapsing the datacube between 2.29 and 2.33^im and over- 
lay this on the HST Jgi4-band image in Fig 6. We extract 
a spectrum from the datacube around the position of the 
submm galaxy and also show this in Fig. 6. The Ha red- 
shift of comp onent A is in excellent agreement with the UV 

redshift from lChapman et alJ ]2005l) . 

The rest-frame UV spectroscopy from lChapman et alJ 

(2005) is measured from absorption lines, yet the high 
[Nii]/Hq emission line flux ratio (0.6+0.2) we see may be 
mor e consistent with either high metallicity or AGN activ- 
ity llVeilleux fc Osterbrockiri987f) . We note that this galaxy 
is not detected in the archival 0.4-10 KeV XMM X-ray image 
of this field to a flux limit of 6xl0 _15 ergs _1 cm -2 (which 
corresponds to Lx <3x 10 44 erg s _1 ). 



4 DISCUSSION 

In this paper, we have spatially mapped the rest-frame op- 
tical line emission in a sample of six high redshift, lumi- 
nous submm galaxies. We also include in our discussion 
two similar galaxies which also have published IFU obser- 
vations (N2 850.4 and S MMJ 14011+0252; ISwinbank et alJ 
120041: iTecza et al-lliooll . The most striking feature in the 
combined sample of eight submm galaxies is the prev- 
elance of galactic-scale multiple components: five of the 
eight appear to comprise two or more galactic-scale exten- 
sions/components on 1-2" (8-16 kpc) scales. Our IFU obser- 
vations also show that velocity structures can be identified 
within individual components (e.g. on ^0.5", 4 kpc scales). 

Such velocity structures are not only obse r ved th rough 
the nebular e missio n: recently iGreve et al] (120051) and 
iTacconi et al] i2006l) have published mm-wavelength ob- 
servations of twelve submm galaxies. Several of these sys- 
tems show double-peaked CO emission line spectra, with 
velocity offsets of up to 800 kms -1 and line widths of up 
to 1000 kms" 1 (FWHM). The individual CO peaks often 
have similar intensities (Ico 1 /Ico 2 =0.7+0.1 with a range 
of 7001/^002=0.3-0.9) suggesting that many of the submm 
galax ies could contain massive gas disks dCenzel et al] 
2003). This interpretation is supported by examples of 
ULIRGS with gas rich disks in the local Universe (e.g. 
NGC6240: ITacconi et ai]ll999l) . However, the obvious in- 
teracting/merger nature of many of these distant submm 
galaxies from our HST and IFU observations strongly sug- 
gests that the velocity offsets we measure arise from merging 
systems. The IFU observations also show that the compo- 
nents have large differences in their energetics and AGN con- 
tributions over 2-20 kpc scales, interpretation is supported 
by recent observations of 

4.1 Mass Estimates 

Constraining the stellar, gas and dynamical masses of 
submm galaxies is crucial if we are to understand their 
evolution. Recent dynamical masses for SMGs from CO 
spectroscopy with millimetre interferometers indicate a dy- 
namical mass of log(Mdj/ n ) = 11.1 ± 0.3 and gas masses 
of log(M oas ) = 1 0.5 ± 0.2 within th e central ~ 10 kpc 
jGreve et alJl2005t ITacconi et al.ll2006l) . I n addition, esti- 
mates of the stellar masses for submm galaxies have been 
given by iBorvs et al] ]2005h . These are derived by fitting 
simple stellar population models to the spectral energy dis- 
tributions of submm galaxies across the restframe ultravi- 
olet, optical and near-infrared wavebands. The inclusion of 
restframe near-infrared data into these fits improves their in- 
sensitivity to dust extinction and thus provides much more 
reliable measures of the luminos ity of the ste ll ar pop ulations 
in these systems. In this way, IBorvs et al.l (I2005T) derived 
a median stellar mass of log(M*) = 11.4 ± 0.4 for a small 
sample of spectroscopically identified submm galaxies in the 
GOODS-N region. 

We are also able to place limits on their dynamical 
masses in this population using the projected velocity offsets 
and spatial separations in our IFU observations. In Fig.|7|we 
construct a velocity offset versus spatial offset diagram for 
all submm galaxies which show multi-components in their 
spatially resolved spectra, either from IFU spectroscopy 
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(this work; ISwinbank et all 120051; iTecza et alJ I2004T): re- 
solved CO spectroscopy jGreve et all 120051: I T^cconi et alJ 
1200ft : or longslit observations ijSwinbank et, al]l2004h . We 
note that no two data-points in Fig. |7| represent the same 
sub-mm galaxy: the only multi-component SMG which has 
both CO and IFU observations (N2 850.4) is represented by 
an IFU data-point in the figure (the measured kinematics 
in N2 850.4 from IFU and CO observations give c onsistent 
results and ar e discu ssed in lSwinbank et all (l2005Tl : see also 
iTacconi et all J2006l) 1. 

To investigate the masses of these systems we assume 
that the velocity offsets between the different components 
eventually transform into motions in a single, relaxed po- 
tential well, i.e. they reflect virial motions. We therefore 
overlay the expected projected velocity dispersion as a func- 
tion ofj^adius_Jorjmel^ NFW pro- 
file iNavarro et al.ll9 97: Lokas & Mamon 2001J). We assume 
a concentration of c = 7 and a virial radius of 200 kpc 
- typical of a massive elliptica l galaxy iMamon fc Lokasl 
120051 : iRomanowskv et alJ l2003t) . Interpreting the velocity 
offsets in this way, we derive a median dynamical mass of 
5±3xl0 n Mq within the virial radius for these galaxies. 

Before we compare the dynamical mass estimates with 
the stellar and gas masses, we note first that there are clearly 
uncertainties in applying the virial theorem and attempting 
to measure the mass of a (non-stationary) merging systems. 
The following caveats must be considered: 

1) Whilst these are merging systems and out of equilib- 
rium, the virial theorem may still be applied if an extra term 
is included to reflect the moment of inertia of the system 
(d 2 I/dt 2 ): such that 2T+W=d 2 I/dt 2 . Since the moment of 
inertia is strongly weighted towards large radii (lex J r 2 dm), 
relatively minor parts of the galaxy (such as tidal tails) can 
have a large impact on / and hence the virial mass estimate. 
However, in many of these systems, the velocity dispersions 
(and star-formation activity) of the individual components 
(and the fact that two SMGs appear to comprise double 
AGN) suggest that we are not measuring velocities for "mi- 
nor" components, but rather these are likely to be the dom- 
inant components of the merger. 

2) The virial theorem is strictly only valid if the com- 
ponents we observe follow the distribution of mass in the 
galaxy. Since the components we measure are the likely nu- 
clei of the progenitor galaxies they will be more concentrated 
(due to dynamical friction) than the overall mass distribu- 
tion and may produce an underestimate of the dynamical 
mass. 

3) The virial theorem should only be applied if the 
merging components have similar mass ratios. The velocity 
dispersions (and star formation activity) of the individual 
components we measure are often comparable arguing that 
the mass ratios of the merging components are similar. It is 
also useful to note that recent observations of local (z <, 0.1) 
ULIRGs which have shown that galaxies with L £ 10 12 Lq 
typically comprise systems two (or more) interacting c om- 
ponents with mass ratios 1:1 to 1:3 l|Dasvra et al.ll20(Ml . 

To gauge the validity of the mass estimate from the 
virial theorem, we can make a further estimate of the dy- 
namical masses: if we take the velocity offsets and assume 
these are point masses (i.e. assume M = Av 2 / ArG), we 
derive masses of 2.5±1.5xl0 M©. We view this estimate 
as more uncertain than that using the virial estimator since 



this estimate requires placing all of the mass of the pro- 
genitors into point masses on circular orbits. The estimate 
based on the point-mass model is within the range we quote, 
and so we consider 5±3xl0 n Mq a reasonable estimate of 
the mass and its uncertainty in the typical submm galaxy 
in our sample. However, we note that these estimates would 
benefit considerably from (high angular resolution) adaptive 
optics IFU observations coupled with a theoretical study of 
the distribution of detectable emission line gas within dusty, 
merging systems at high redshifts. 

4-1.1 Comparison with Stellar Mass Estimates 

iBorvs et alJ <l2005h estimated restframe 2.2-/im luminosities 
for a sample of spectroscopically-confirmed submillimetre 
galaxies in the GOODS-N region from broad-band pho- 
tometry covering UBVRIz JK+IRAC 3.6/4.5/5.8/8/im - 
which spans the restframe UV, optical and near-infrared 
at the redshifts of their sample, z = 0.56-2.91. They find 
Mk = —26.8 ± 0.4 from their sample of 10 galaxies at 
z > 1.5. They then adopted a typical light-to-mass ratio for 
the stellar population of La-/M=3.2, based on an estimated 
mean age of ~ 250 Myrs from fitting simple stellar popula- 
tion models to the broad-band SEDs, a nd so derive typica l 
stellar masses of log(M„) = 11.4 ± 0.4 ijBorvs et al.ll2005ft . 
This is a significan t stellar mass for a galaxy at z ~ 2 and so 
IBorvs et alJ i2005l) discuss in detail the potential uncertain- 
ties in their estimates, in particular whether it is possible to 
reduce the stellar masses. The most obvious way to achieve 
this would be to adopt a higher Lx/M ratio, equivalent to 
assuming a younger age for the stellar population. However, 
to achieve an order-of-magnitude reduction, i.e. a light-to- 
mass ratio of Lx/M~ 30, wou ld require a typical age of 
~ 10 Myrs l|Leitherer et al]ll999T ). Given the observed space 
density of submillimetre galaxies, such very young ages and 
hence high-duty cycles imply very large volume densities for 
the progenitor population (~ 10 -3 Mpc~ 3 ) which is at least 
an order of magnitude greater than the estimates of the 
masses and clustering of this population ijBlain et, alJl2004l : 
IChanma.n et aIll20oA. 

It is also possible that AGN contribute to the measured 
2.2-/im luminositi es, leading to an o ver-estimation of the 
stellar masses. As lBorvs et alJ l|2005h point out, the shape 
of the spectral energy distributions for typical submillimetre 
galaxies is well-described by a stellar model and indeed they 
remove from their analysis the small number of galaxies for 
which an AGN is potentially dominant in the restframe near- 
infrared. Nevertheless, we can take an extreme view point 
and ask if AGN produce all of the luminosity in the red- 
dest bands in these sources, how much would that influence 
the estimated stellar luminosity? We use the average near- 
infrare d spectral index for QSOs from lSimpson fc Rawlinesl 
( 2000), a = 0.98, assume that all of the observed 8.0-fim lu- 
minosity is arising from the AGN and estimate its contribu- 
tion to the luminosity at a restframe wavelength of ~ l.Qfim 
(where the stellar contribution should peak). Using this set 
of extreme assumptions, we find a median reduction in the 
stellar luminosity of only 4 x , suggesting that a more realis- 
tic estimate taking into account the stellar contribution at 
8.0/im, is likely to yield a ch a nge o f ^ 2x, within the error 
range quoted by|Borxa_et_ajJ <l2QQ5l) . It is also worth noting 
at this point that TAlexancteretHilT i2005h have shown that 
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Figure 7. Velocity offset versus spatial offset for the multiple 
components within submm galaxies from IFU, longslit or CO ob- 
servations. The dotted lines show the expected profiles as a func- 
tion of halo mass for a galaxy with an NFW profile and a con- 
centration parameter typical of elliptical galaxies. Treating the 
observations as an ensemble, we derive a dynamical mass for the 
SMGs of M=5±3xlO n M (shown by the solid line). For com- 
parison, the dashed line shows the profile for the minimum stellar 
mass derived from Borys et al. (2005). 



while AGN are ubiquitous in SMGs, their typical contribu- 
tion to the bolometric emission is & 10% suggesting that 
any contamination to the 2.2/xm flux from AGN activity is 
minimal and unlikely to bias our result. 

Finally, we note that recent stellar population synthesis 
modelling by Maraston et al. (2006) has shown that ther- 
mally pulsating asymptotic giant branch (TP-AGB) stars 
in post starburst galaxies can contribute significantly to the 
rest-frame K-band luminosity on £ 200-1000Myr timescales. 
Indeed, the pre dicted stellar light-to-mass ratios at 250 Myr 
is Lk/M~ 6-8 iMarastonlll998l : iMaraston et alJl2006l) and 
therefore may act to reduce our estimated stellar mass by 
a factor of up to 3 x if the luminosity weighted stellar pop- 
ulations in the SMGs are dominated by stars with ages of 
£ 200-1000Myr. However, in the case of an on-going star- 
burst (as is relevant for SMGs) the importance of TP-AGB 
phase is likely to be lessened and we estimate a maximum of 
L K /M ~ 5 at 450 Myr and more typically h K /M g 4 for the 
likely ages of SMGs, suggesting that this phase is unlikely 
to bias our results significantly. 

Hence, unless there is a conspiracy between a number of 
factors, it appears unlikely that the estimated stellar masses 
for the submillimetre population are over-estimated by more 
than a factor of 2-3 x. We conclude therefore that the typi- 
cal stellar and dynamical masses for submillimetre galaxies 
are likely to be similar in magnitude. We estimate a mass 
ratio of Md y n,/M* ~ 2, with a considerable uncertainty. This 
suggests that, as with local luminous ellipticals, the central 
regions of these galaxies are baryon-dominated. 



4.2 Space Densities and Comparison with Local 
Populations 

In order to test the connection between SMGs and local 
galaxy populations we can evolve their stellar populations 
to the present day and determine whether their descendents 
would be consistent with scaling relation s for local g alax- 
ies, s uch as the Faber-Jackson relation jFaber fc Jacksonl 
1976). This will test claims that SMGs evolve into mas - 
sive ellipticals at the present day (e.g. iLillv et all ll999). 
The velocity offsets we measure in H4. II implies an average 
velocity dispersion of approximately a — 240 ± 80kms _1 . 
The Faber-Jackson (FJ) relation for local early-type galax- 
ies indicates that this dispersion c orresponds to an r-b and 
magnitude of M r = -22.7 ± 0.5 llBernardi et al.ll2005l) or 
Mk = —25.3 ± 0.5 (assuming (r — K) = 2.6, typical for 
an elliptical galaxy at the present day). Is the expected 
fading of the stellar populations within the SMGs, with 
Mk = —26.8 ± 0.4 at z ~ 2, sufficient to reconcile them 
with the z = FJ relation? 

The key to estimating the degree of fading of the SMGs 
is to determine the likely ages of the stellar populations. 
Since the starburst lifetimes in SMGs depends on the mass 
of the molecular g as reservoir and the star-formation rate, 
iGreve et all J2005I) used CO (3-2) transitions for 18 SMGs 
(and assuming a CO-H2 conversion factor) to infer gas 
masses, obtaining M(H2) ~ 3 x 10 10 Mq. Recently, Hain- 
line et al. (2006) have reported CO (1-0) observations of an 
SMG using the Green Bank Telescope and derive a H2 gas 
mass four times larger than that from the CO (4-3) transi- 
tion of the same source (implying that J upP er >3 transitions 
of CO significantly underestimate the mass of cold, diffuse 
molecular gas in SMGs). This suggests that the typical gas 
masses of SMGs are likely to be M(H 2 ) ~ 1 x 10 11 M . Star 
formation rates can be derived from the bolometric emission 
of the source, assu ming this is powered by star formation. 
IGreve et all <200SD used 850 /im and 1.4 GHz fluxes and as- 
sumed the far-infrared-radio correlation to determine bolo- 
metric luminosities for their sample. They then assumed a 
conservative contribution of 50% to the bolometric luminos- 
ity from a AGN, and so estimated a median star-formation 
rates of 7OOM0yr~ 1 for the galaxies in their sample (inte- 
grating the Salpeter IMF between 1 and 12OM yr _1 ). Two 
recent results influence this conclusion, b ut ess entially leave 
it unchanged. Firstly, lAlexander et alJ i2005l) have shown 
that while AGN are ubiquitous in SMGs, their typical con- 
tribution to the bolometric emission is & 10%, which would 
increa se the star-formation rate estimated by IGreve et alJ 
(2005) by ~ 2x. However, Kovacs et al. (2006) have used 
350 fim photometry to improve the estimates of the bolo- 
metric luminosities from the far-infrared SEDs of SMGs and 
show that the star-formation rates may be overestimated by 
a factor <; 2x, cancelling this increase. Hence, our current 
best-estimate of the gas consumption timescales in SMGs 
is: tsmg ~ M{H 2 )/SFR ~ 1 x 10 11 M /700 Meyr" 1 -^ 
150 Myr, although this estimate clearly has significant un- 
certainties. The reader should note that a burst with this 
duration and star-formation rate can account for all of the 
observed stellar mass in a typical SMG, suggesting that the 
activity we are witnessing is likely to be the origin of the 
majority of the stellar mass in these galaxies iSmail et al.l 
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We can next compare these burst lifetimes with those 
required in order for the SMGs to lie on the FJ relation at 
the present day using stellar synthesis modelling. For SMGs 
with z > 1.5, th e aver age rest-frame if-band magnitude 
from iBorvs et alJ J2005I) is Mk = -26.8 ± 0.4. We start 
by calculating the fading of a simple inst antaneous burst 
to th e present day. We us e starburst99 llLeitherer et alJ 
Il999l) with a Salpeter IMF JSalpeterlll955h and an instanta- 
neous star-burst which predicts that between 100 Myr and 
10 Gyr there should be approximately AMk = 3.4 mags of 
fading suggesting the present day descendants should have 
Mk ~ —23.4 ± 0.4. This is approximately 2 magnitudes 
fainter than the FJ relation, given the expected velocity dis- 
persion as found in H4.ll Similarly, an instantaneous burst 
using the models of Maraston et al. (2006) suggests that 
the fading between 100 Myr and 10 Gyr is approximately 
AMk = 2.1-2.4 (depending upon the metallicity) - which 
is still £ 0.6 magnitudes fainter than the FJ prediction, given 
our expected velocity dispersion. 

However, these instantaneous burst models are unphys- 
ical for SMGs which are still forming stars at the current 
epoch. So, we derive the fading from a more realistic model: 
a constant burst of star-for mation of duration 1 00, 200 
and 300Myrs using PEGASE (|Le Borgne et al.ll2004h . For a 
100 Myr burst (observed at 50 Myr, since we expect on av- 
erage to see SMGs half way through their active phase) the 
fading between the epoch of observation and the present day 
is AMk = 3.2, ~ 1.5 magnitudes too faint compared to the 
FJ relation for the velocity dispersion. Taking longer burst 
durations of 200 and 300 Myrs (observed at 100 and 150 Myr 
respectively) the fading is AMk = 2.7 and 2.3 magnitudes 
respectively - approximately 1.2 and 0.8 magnitudes off the 
FJ relation. Thus for our preferred burst lifetime, 300 Myrs, 
observed half way through the burst, the predicted present- 
day magnitude is Mk ~ —24.5 ± 0.4 compared to the 2 = 
value of Mk = —25.3 ± 0.5 from the FJ relation. 

These predictions are of course sensitive to the choice of 
IMF, changing from a Salpeter to a Kroupa IMF and for a 
300-Myr burst the predicted fading reduces to AMk = 1.7, 
indicating a present-day luminosity of Mk ~ — 25.1 ± 0.4 
and placing the SMGs within ~ 0.2 magnitudes of the 2 = 
FJ relation. We note that such a change would be at odds 
with the "top-heavy" IMF needed in semi-analytic models of 
galax y formation to rep roduce t he basic properties of SMGs 
(e.g. feaugh et al]l2005T: see also lBlain et al.l|l999blah . 

Hence, whilst there are a number of uncertainties in es- 
timating both the starburst lifetimes and the fading of SMGs 
to the present day, longer burst durations (~ 300 Myrs) are 
favoured in order for SMGs to evolve towards the scaling 
relations for elliptical galaxies at the present day. Burst 
timescales of this length are also consistent with the ex- 
pected time to consume the gas reservoirs in these galaxies 
and with the build up of the present stellar mass in the 
systems. 

We can also investigate how the SMG burst lifetimes 
influence the space densities of their expected descendants 
at the present day. Th e full SMG redshift distribution from 
man et alJ 1)20051) is well fit by a Gaussian with a mean 
redshift of < 2 >= 2.2 and a with, a z ~ 1.3 (accounting for 
incompleteness in the spectroscopic desert). Thus the SMGs 
populate an epoch of about T b s = 1.5 Gyr in duration. Since 
"^psMGTobs/Ttmrst (where <f) is the comoving space density 



of the descendant population, Tb urs t is the burst lifetime, 
and psmg is the observed space density) an estimate of 
Tburst allows us to infer <j). We first estimate the observed 
space density (psmg) by t aking the observed surf ace den- 
sity of >5mJy SMGs from IChapman et al.1 J2005I) (~1000 
per square degree) which suggests that between 2=0.9-3.5 
the volume density should be ~ 2.5 x 10 _5 Mpc -3 (this cal- 
culation assumes that 67% of the SMGs lie within the la 
limits and includes the incompleteness correction in both 
the redshifts and radio-detection limits). Using our esti- 
mated lifetime of Tb urs t ~ 300 Myr, the corrected space den- 
sity (assuming a galaxy only goes through one SMG phase 
in its lifetime) is therefore 2.5 x 10~ 5 x 1.5 Gyr/300 Myr ~ 
1.3 x 10 -4 Mpc -3 for a population with predicted present 
day isT-band absolute magnitude of Mk ~ —25.1. Integrat- 
ing dow n the 2 = luminosit y function of early-type galax- 
ies from ICroton et al.l (|2005|) (see also Wake et al. 2006), 
to a median magnitude of Mk ~ —25.1 we determine a 
space density of 9 x 10 -5 Mpc -3 (we note that using the 
LRG luminosity function from Wake et al. 2006 we ob- 
tain 1.1. x 10~ 4 Mpc~ 3 ). Thus if SMGs have burst dura- 
tions of approximately 300 Myr, as indicated by their gas 
consumption timescales, then their predicted fading would 
place them on (or near-to) the present day FJ relation and 
the observed space density of this population at 2 ~ 2 (cor- 
rected for the duty-cycle of the bursts) would account for 
the formation of the whole population of early type galaxies 
with Mk < -25.1, 3L* K ) seen at 2 ~ 0. 

Further evidence for the evol ution of SMGs into ma ssive 
ellipticals has also be presented in lTacconi et alJ <|2006T) who 
used high resolution millimeter imaging to derive matter vol- 
ume and surface densities of ~ 100 cm -3 and 5000 Mqpc -2 
for the SMGs, comparable t o those in massive local ellipti- 
cals and Sa bulges. Similarly. iBlain et al.l i2004l) provide ten- 
tative evidence for clustering of SMGs at a level consistent 
with the progenitors of > L* elliptical galaxies. These pro- 
vide additional support for our conclusion that SMGs can 
be connected through a self-consistent evolutionary model 
to the local luminous galaxy population, which reproduces 
both the mass-luminosity relationship and luminosity dis- 
tributions of local elliptical galaxies. 



5 CONCLUSIONS 

In this paper we have studied the rest-frame optical emission 
line structures and dynamics of six high redshift, luminous 
submm galaxies. Including two sources from the literature 
with similar observations, we find that of the eight submil- 
limetre galaxies with resolved spectroscopy, at least five ap- 
pear to comprise two or more dynamical sub-components. 
The average velocity offsets between these components is 
~ 180 km s -1 across a projected spatial scale of 8kpc. The 
obvious merging/interacting nature of these systems sug- 
gests there are analogous to the multi-component nature of 
the, typically less bolometrically luminous and slightly more 
compact, ULIRGs in the local Universe. 

Since our IFU observations allow us to disentangle 
AGN and star-burst-like components (e.g. from [Nil] /Ha 
flux ratios), it is perhaps significant that two of the eight 
submm galaxies shows possible signs of AGN activity in two 
spatially-resolved components on 1-2" (10-20 kpc) scales. 
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Deep X-ray surveys of the Chandra Deep F i eld Nort h have 
produced similar results: lAlexander et alJ i2003al) report 
that two of seven submm galaxies in this region are indi- 
vidually associated with pairs of X-ray sources on similar 
spatial scales (approximately one galactic diameter). Since 
the chance of a false association is g 1%, it may be that we 
are witnessing the interaction or merging of AGNs in these 
sources (a low redshift example of t his binary AGN activ - 
ity is seen in the ULIRG NGG 6240: iKomossa et"ai1l2003h . 
The double AGN in some of these systems, suggests merg- 
ers may be responsible for disturbing the gas orbits in both 
components, hence fuelling subsequent AGN growth which 
curren tly lags behind a significant pre-existing st ellar popu- 
lation iAlexander et a l. 2003a; Ismail et al]l2003f) . 

By combining the spatially-resolved kinematic informa- 
tion from IFU, CO and longslit observations and assuming 
the velocity offsets we measure are eventually transformed 
into velocity dispersions in a pressure-supported galaxy, we 
derive dynamical masses of 5±3xlO n M for typical lu- 
minous submillimetre galaxies. This is si milar to the dy- 
nami cal masses found using Ha line widths 1^ winbank et alJ 
2004) an d resolved CO spectroscopy iGieve et ajj 2005; 
iTacconi et alll2006h and a factor of ~ 2x grea t er tha n the 
recent stellar mass estimates from lBorvs et all ((2005) , sug- 
gesting that the central regions of these galaxies are baryon 
dominated. 

We combine our dynamical measurement with esti- 
mates of th e stell ar luminosities for this population from 
iBorvs et alJ i2005f) . We use these to constrain their likely 
evolution and so compare the properties expected for the 
descendents of submillimetre galaxies with z — galaxy 
populations . We update the g as depletion time scales es- 
timated bv lGreve et al.l i2005l) using recent results on the 
gas masses, bolometric luminosities and AGN contributions 
for this population, deriving typical gas depletion times of 
~ 150 Myr. Adopting a typical burst duration of twice this 
age, ~ 300 Myr, we predict the evolution of their stellar pop- 
ulations and compare these to the Faber- Jackson relation for 
local ellipticals. For a ~ 300 Myr burst with a Kroupa IMF, 
we find that the expected fading would result in the descen- 
dents of the SMGs lying within ~ 0.2 mags of the FJ relation 
at their measured velocity dispersion. The same estimate of 
the lifetime can be used to correct the apparent space den- 
sity of SMGs (assuming a galaxy only goes through one SMG 
phase in its lifetime) to derive a space density for the descen- 
dents of this population of ~ 1.3 x 10 _4 Mpc -3 for a popula- 
tion with z — absolute luminosities of Mk ~ —25.1. Hence 
if SMGs have burst durations of ~ 300 Myr, as indicated 
by gas consumption timescales, then their predicted fading 
places them close to the present-day FJ relation. Moreover, 
the observed population at z ~ 2 is very comparable to that 
of luminous ellipticals with Mk S —25.1 (~ 3L* K ) at the 
present day. 

Our observations show the power of using integral field 
spectroscopy around the rest-frame optical emission lines 
to probe the dynamics and power sources of these galaxies 
on ^5 kpc scales. With forthcoming adaptive optics integral 
field spectroscopy on 8- and 10-m telescopes (e.g. OSIRIS 
on Keck, NIFS on Gemini and SINFONI on the VLT), we 
will soon be able to probe the dynamics, metallicities and 
star-formation properties on &0.1" (sub-kpc at z = 2) scales. 
Furthermore, the increased light grasp of large aperture tele- 



scopes will allow us to spatially resolve and map the ab- 
sorption line spectra of these galaxies (including important 
diagnostic lines such as Can H&K, CH G-band, H<5 and 
MgB), allowing us to accurately probe the stellar contents 
and metallicity gradients, and so investigate the stellar ages 
and populations produced in these extreme systems. Such 
observations will allow us to probe further the timescales 
involved in the encounters as well as how and when these 
galaxies form most of their stars and thus the precise relation 
to similarly massive galaxies at the present day: luminous 
ellipticals. 
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